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■BSTRACT 


We  consider  ^ special  class  of  charac- 
teristic first-order  initial  value  problems 

P(x  ,X  ,^,P1#P2)  = o .S> The  initial  value 

problem  arises  in  the  asymptotic  solution 
of  parabolic  and  elliptic  equations.  The 
problem  is  characterized  by  a singular, 
characterisitc  initial  manifold.  Namely, 
initial  data  is  given  on  a characteristic 
curve.  The  characteristic  curve  is  also 
sinqular  in  that  there  is  a point  on  the 

/'"2  2 

initial  manifold  where  ( F + F \=  0 . 


■We-prov^  £hat  such  problems  have  unique 
solutions^*1  The  theorem  also  has  an 
interesting  interpretation  in  terms  of 
the  calculus  of  variations. 


Section 

Section 


SECTION  1 
INTRODUCTION 

In  this  note,  we  will  discuss  a characteristic  initial 
value  problems  (IVP)  that  is  more  singular  than  the  usual  IVP, 
but  has  a unique  solution.  The  problems  arise  in  the  asymptotic 
solution  of  elliptic  and  parabolic  equations  connected  with 
stochastic  dynamical  systems  (Ludwig,  1975;  Mangel,  1977;  Mangel 
and  Ludwig,  1977) . An  example  of  the  type  of  initial  value  problem 
of  interest  is: 


b1i^i  - -j-  (i|/)  = 0 


t|i  = vl)  on  S ; f (^  ) = 0 

r o o 


S parameterized  by  x1  = b1 


(1) 


(2a) 


(2b) 


In  equation  (1) , subscripts  indicate  partial  derivatives  and  repeated 
indices  are  summed  from  1 to  n . The  problem  posed  by  (1)  - (2) 
is  a characteristic  initial  value  problem,  since  S is  a charac- 
teristic curve  of  (1)  (Courant,  1962) . In  the  stochastic  problems, 

S is  also  a singular  manifold.  Namely,  there  exits  one  point  Q€S, 


such  that 


b (Q)  = 0 for  all  i 


(2c) 


i 


Under  the  additional  condition  (2c) , it  was  shown  (Mangel,  1977) , 
that  (1)  - (2)  has  a unique  solution.  In  this  note,  we  will 
generalize  the  above  result.  An  analogy  to  the  partial  differential 
equation  is  presented  in  section  2,  where  we  discuss  an  ordinary 
differential  equation  with  a singularity.  In  section  3,  the  main 
theorem  is  given.  In  the  proof  of  the  theorem,  we  use  the  singularity 
of  the  initial  manifold  to  construct  the  unique  normal  derivative 
of  <Mx)  on  S . Once  the  normal  derivative  is  known  uniquely, 
it  is  possible  to  calculate  i|>(x)  off  the  initial  manifold  using 
the  method  of  characteristics.  A comparison  with  the  theory  of  the 
standard  initial  value  problem  is  given.  In  section  4,  we  discuss 
the  theorem  of  section  3 in  the  setting  of  the  calculus  of  varia- 
tions. Finally,  in  section  5,  we  return  to  the  example  given  by 
equations  (1)  - (2a, b) , and  show  how  the  unique  solution  can  be 
obtained  for  the  case  f(i|>)=^. 


SECTION  2 

ORDINARY  DIFFERENTIAL  EQUATIONS  WITH  A SINGULARITY 


A simple  analogy  to  the  problem  presented  in  section  1 is 
the  following  ordinary  differential  equation: 


b(x)  ||  + f (x,  ty)  = 0 


We  assume  that  there  exists  exactly  one  point  xp  such  that 


b(x  ) = 0 


b (xp)  ? 0 


(3) 


(4) 

(4') 


We  require  that  ^ is  regular  at  Xp  . Since  b(Xp)  = 0 , for  \p  to 
be  regular: 

f (x  » ip  (x  ) ) = 0 (5) 

r r 


We  assume  that  (5)  defines  <Mxp)  uniquely.  We  note  that  the 
regularity  condition  imposes  initial  data  on  the  equation.  The 
derivative  d /dx  at  xp  is  obtained  by  differentiating  (3) : 


r 


solution,  equation  (3)  has  the  unique,  regular  solution 


^ (x) 


-I 


f(t,  ip(t))dt 
b (t) 


+ ^(xp) 


(7) 


We  obtained  a unique  solution  without  specifying  any  data  for  ; 
the  only  requirement  was  that  it  be  regular.  The  singularity 
imposes  data  on  the  equation.  Many  examples  can  also  be  found  in  the 
theory  of  second  order  differential  equations  (e.g.,  Bessel's  equation). 

We  will  show  that  a similar  phenomenon  occurs  with  the  problem 
(1,  2).  Namely,  if  we  require  a regular  solution,  then  the  singularity 
on  the  initial  manifold  will  force  a unique  value  for  the  normal 
derivative  of  i/>  (corresponding  to  di^/dx  evaluated  at  x ) . 

It 


r 


SECTION  3 

UNIQUE  SOLUTIONS  OF  CHARACTERISTIC  INITIAL  VALUE  PROBLEMS 


In  this  section,  we  prove  the  following  theorem.  Let 


1 2 


F (x  ,x  =0  be  one  first  order  partial  differential 

equation  . Let  Cauchy  data  ij;  = ^(t)  on  a manifold  S ^ be  given. 
The  initial  surface  S is  parametrized  by 


(8) 


dx 


dt 


= b1 (x) 


(9) 


We  assume  that: 

i)  S is  a characteristic  manifold,  i.e., 


A = 


F F 

'P. 


= 0 


on  S; 


(10) 


ii)  The  mainfold  S has  exactly  one  singularity,  i.e.,  there 
exists  one  point  QgS  such  that 

F^1(Q'  ^l'  ’^2)  + F\p2(Q'  *1'  ^2)  = 0 

for  all  values  of  \l>2  • 

iii)  Let  n denote  distance  normal  to  S . We  assume  that 

«>  v ■ 0(n“> 


(id 


n-0 


(12) 


I 


r 


b)  U”  F*  *1  % *2  ' 0(n  > 

n-0  yl  r2 


c)  Let  x denote  distance  normal  to  S.  Then 


(12) 


lim  b1(x)  = 0 (na) 
n-*0 


Under  these  assumptions,  the  initial  value  problem  has  a unique 
solution.  Namely,  it  is  possible  to:  i)  calculate  a unique  value 
of  the  normal  derivative  of  ip  on  S ; ii)  integrate  the  charac- 
teristic equations. 

PROOF 

We  prove  parts  i,  ii  of  the  conclusion  separately. 

i)  Calculation  of  \b  on  S 
_n 

As  is  usually  done,  we  assume  that  \p  has  as  many  smooth 

derivatives  as  are  needed.  Without  loss  of  generality,  we  assume 

1 2 
that  x represents  distance  normal  to  S and  x represents 

distance  along  S . We  set  the  origin  at  the  singularity  Q. 

The  tangential  derivative  of  \p  can  be  obtained  by  differentiation 

of  the  initial  data  i^(t)  . 

We  differentiate  (8)  with  respect  to  x1  and  obtain 

Fxl  + Vl  + + = 0 

At  the  singularity  Q , equation  (11)  holds,  so  that  (13)  becomes 


FX1(Q,  ip,  iplt  <l>2)  + F (Q,  \p,  (|>ir  if^)^  = 0 


(13) 


(14) 


l 


Li 


In  equation  (14),  Q,  41,  and  ip2  (the  tangential  derivative)  are 
known.  We  assume  that  (14)  yields  a unique  value  for  ^(Q)  . 
Next,  we  use  the  characteristic  condition  (10) , which  can  be 
rewritten  as 


= F, 


= F., 


We  use  (15)  to  express  F in  terms  of  F 


V + Vi + -j§  (bl*n  + ^i’  - 0 

1 > •-  t \ 1 1 1 *\  / ^ .i. 


Equation  (13)  becomes 


In  light  of  (9),  d/dt  = b 9/3x  , so  that  (16)  can  be  rewritten  as 


^2-  {F  1 + F.i|>,  ) 
F,  1 x ^rlJ 


The  characteristic  condition  (10)  indicates  that 


(17a) 


Hence,  by  assumption  iii)  b /F.  is  bounded  away  from  zero  and 
finite.  Consequently,  ^ will  be  a regular  function  on  S , given 


'i(t>  -l 


2 

2-iSI  (FX1  + f *j)ds  + t1(a) 
*2 


§ 


We  assume,  as  is  usually  done,  that  (18)  is  a meaningful  solution 
of  (17).  We  note  that  if  the  singularity  were  not  present,  then 
tji^(t)  would  not  be  determined  uniquely.  Instead,  we  would  have 
obtained  a one-parameter  family  for  i|^  (see  remark  2,  also). 
ii)  Integration  of  the  Characteristic  Equations 
The  characteristic  (or  ray)  equations  corresponding  to  (8) 


dx_  = F 
dx 


= F, 


-a?  *i  + 


dx  , 
TF  *2 


f(x,  ^) 


- (p,  F.  + F , ) 
x 


k = 1,  2 


In  the  last  equation  in  (19) , we  have  used  the  notation  that 
= p^  . Since  the  initial  manifold  S is  characteristic,  we 
have  that 


= 0 


Equation  (20)  indicates  that  it  is  not  possible  to  integrate 
(19)  to  move  off  of  the  initial  manifold.  In  order  to  use  the 
characteristic  equations,  we  construct  a new  manifold  , which 

is  not  characteristic.  The  new  manifold  is  characterized  by 
^ = tj^(t)  + y on  S^,  . We  construct  by  a Taylor  expansion: 


i 


{(x1,  x2):  x1  = y/^fx2)  + 0(y2)} 


We  now  consider  the  new  initial  value  problem:  solve  (8)  with 


(21) 


initial  data  given  on 
ray  equations  (19),  on  S 


We  need  to  give  initial  data  for  the 
We  have  shown  that  ip(0;y)  = \p^  + y , 


where  the  first  argument  corresponds  to  the  ray  parameter  t and 
the  y denotes  that  we  are  on  S 


Also,  x^OjY)  = y/^i  + 0(y“) 


Clearly,  x2(0,y)  = x2  and  ^(OfY)  = 0(y)  + i|>2  • 

Finally,  we  have  ^^(0;y)  = + 0(y)  , where  is  the  value  of 

^ on  S. 


We  obtain 


dx 


dx 


x=0 


= F (■— 


Yf  ^ ^2)  + 0 (Y  +a) 


dtp 

dx 


= f ( x , y)  + 0(y2+“) 


(22) 


(23) 


<7=0 


We  now  reparametrize  the  ray  equations  by  introducing  o 
where  a is  given  in  (12) . Hence,  we  obtain 

F 


a 

TY  , 


dx' 


da 


<Pi 


a= 0 


diji 

da 


+ o (y) 


= LlXj-ll  + o (y) 


(24) 


(25) 


a=0 


We  note  that  in  some  cases,  equations  (22,  23)  may  have  to  be  re- 
parametrized separately.  We  now  let  y-»0  , so  that  S collapses 


9 


into  the  original  manifold  S . We  obtain 


lim 

dx1 

lim 

F.  1 

Jl_  * 

Y+0 

da 

Y+0 

a r 

a=0 

Y 

lim 

dip 

lim 

f (x,y) 

Y^O 

da 

rr-r\ 

Y-+0 

a 

Y 

(26) 

(27) 


Consequently,  we  can  calculate  x^  and  ip  off  the  initial  manifold, 

2 

by  using  (26,27)  and  the  equations  for  x and  from  (19). 

Thus,  in  a neighborhood  of  S , <J;  is  uniquely  determined. 


REMARKS 

1.  An  intuitive  description  of  the  reparametrization  is  as 

follows.  For  each  y^O,  equations  (19)  are  solved,  so  that  we 
obtain  rays  emanating  from  . We  denote  this  set  of  rays  by 

{R^j  . As  y-*0,  equations  (26,27)  indicate  that  {R^}  converges 

to  {Rq},  raYs  that  appear  to  emanate  from  the  initial  manifold  S . 

2.  In  the  usual  treatment  of  the  Cauchy  problem  (Courant,  1962) , 
it  is  assumed  that  on  S 


0 . 


(28) 


Hence,  the  usual  theory  is  not  applicable  to  the  problem  discussed 
in  this  theorem.  The  standard  theory,  furthermore,  indicates  non- 
uniqueness for  characteristic  initial  value  problems.  The  theorem 
given  here  has  extended  the  theory  in  Courant  (1962).  It  is  noteworthy 
that  the  more  singular  problem  has  a unique  solution. 


SECTION  4 


THE  VARIATIONAL  SETTING  AND  HAMILTON- JACOBI  THEORY 


We  now  consider  a first  order  partial  differential  equation 


that  is  independent  of 


For  example,  the  substitution  <(> 


is  a Hamiltonian 


and  the  variational  problem:  choose  <Mx)  such  that 


* ■ CU,X„,X)  1 !.(*<«)).*(»)  Ids 


The  problem  is  solved  when  we  can  construct  level  curves  of  the 
function  <Mx)  (Rund,  1966).  These  curves  are  determined  by  solv- 
ing the  Euler  or  Hamilton- Jacobi  equations: 


In  the  usual  formulation  of  such  variational  problems  (Rund, 1966) 
it  is  assumed  that  for  all  x: 


* 0 


i 1,  . . . , n 


The  theorem  of  section  3 has  as  an  application  the  following 
singular  variational  problem.  Suppose  that 


3p7  (xQ,p)  = 0,  for  all  i 

and  that  9H/9p^  vanishes  nowhere  else.  Also  suppose  that  we 
are  given  exactly  one  extremal  through  xQ  and  the  value  of  <}>  on 
this  extremal.  The  theorem  of  section  3 indicates  that  the 
variational  problem  will  have  a unique  solution. 


I 

i 


A certain  trade-off  in  information  is  involved  in  this 


application.  Namely,  instead  of  (36),  we  require  the  knowledge 
of  one  extremal  curve.  Such  problems  have  not  received  much 
attention  in  the  calculus  of  variations.  However,  exactly  this 
sort  of  problem  arises  in  the  asymptotic  solution  of  certain 
diffusion  equations  (Ludwig,  1975),  (Mangel,  1977),  (Ventcel  and 
Freidlin,  1970). 


SECTION  5 


AN  APPLICATION 


We  conclude  by  showing  how  the  theorem  of  section  3 applies 
to  a special  case  of  the  problem  posed  by  (1)  and  (2) . We  set 
f(^)  = The  initial  data  is  f = 0 on  S.  We  assume  that  Q is 
the  singularity  on  S.  Let  B = (b  * 1 j ) I q > we  assume  that  B 
has  one  real  negative  (X-)  and  one  real  positive  (A+)  eigen- 
value. The  manifold  S is  calculated  by  moving  away  from  Q in 
the  direction  of  the  eigenvector  of  the  negative  eigenvalue  and 
integrating 


dx'*" 

dt 


-b1 (x) 


i = 1,  2 


The  problem  as  posed  is  characteristic.  When  (1)  is  differenti- 

V , 

a ted  with  aspect  to  x and  evaluated  on  S we  obtain 


+ bl*ik 


all 


0 


or 


d*k 

^ 4- 


+ b 


Jt 


0 


Equation  (40)  can  be  used  to  derive  an  equation  for  4*  (since 

n 

= 0 on  S)  . We  find  that  ip  satisfies 


dt**n  . a ,3 

at"  + D*n  - 2 *n 


= 0 on  S 


A A 


The  functions  b,  a are  expressed  in  terms  of  the  original  b1, 

(a1^),  (Mangel,  1977).  In  this  case,  Q corresponds  to  t = 00 

dip 

At  Q,  = 0 for  all  i,  so  that  (41)  yields 


or 


b*n(0°)  -|4'n(") 


ip  (00) 

n 


= 0 


' A 

2b  (°°) 



a (°°) 


Equation  (41)  is  a form  of  Abel's  equation  and  has  the  solution 

s -V2 

*n(t)  = ( f a (S)exp  -?J  b(S')dS'j 


-2.J  b(S')dS'J  dsj 


We  have  constructed  the  unique  normal  derivative  of  on  S. 

We  now  consider  the  ray  equations.  We  switch  to  normal  (N) 
and  tangential  (Y)  coordinates  (on  S ) , with  the  origin  at  Q. 

In  the  vicinity  of  Q,  the  ray  equations  are 


57  = X+N  - aNNpN*  + 0(N2  + Y2) 


15 


(44) 


l 


r 


dY 

dx 


A.Y  - + 0(N"  + Y*) 


dx 

dT 


5 aNNpN^  + 0(N2  + Y2) 


dp 


k 

dx 


b\pi  + ~ 2— Pipj^  + alDpipjpk 


We  introduce  a new  manifold  S ^ , defined  by  ip 


parametrized  by 


S*  = { (N,  Y):  N = y/<I»n(Y)  + o(y2)} 


On  S , equations  (44,46)  become 


dN 

df 


A+y  NN  0 , ,2 

— - a p.„Y  + 0 


T = 0 


N 


PNY  + O(Y-) 


d 

dr 


1 NN  0 ,2. 

" 2 PNy  °(y  ) 


T = 0 


We  reparametrize  by  a ? ry  and  find 


lim  dN 
y->0  cTo 


o=0 


^-aNNp°  + 0 


lim  dx 


1 NN  0 , 


Hence,  the  ray  equations  (51,  52,  45,  47)  can  be  integrated  to 
determine  <Mx)  off  of  S . 

Other  examples  of  such  characteristic  initial  value 
problems  are  given  in  Mangel  (1977). 
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